On a relation between the Bach equation and 
the equation of geometrodynamics 

M.V. Gorbatenko, A.V. Pushkin, H.-J. Schmidt 

Abstract 

The Bach equation and the equation of geometrodynamics are 
based on two quite different physical motivations, but in both ap- 
proaches, the conformal properties of gravitation plays the key role. 
In this paper we present an analysis of the relation between these two 
equations and show that the solutions of the equation of geometro- 
dynamics are of a more general nature. We show the following non- 
trivial result: there exists a conformally invariant Lagrangian, whose 
field equation generalizes the Bach equation and has as solutions those 
Ricci tensors which are solutions to the equation of geometrodynam- 
ics. 
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1 Introduction 



Conformal properties for theories of gravity appear on several different levels. 
The best known example is the conformally coupled scalar field, where the 
conformal property is due to the fact, that in four dimensions, the operator 
□ — R/6 is conformally invariant. 

Two other examples are: First, the Bach equation B a p = 0, (see e.g. 

and and the references there for details, ||, |4j] for further recent 
approaches to the Bach equation, and |5j for an overview on fourth-order 
gravity in the time period 1918 - 1990) whose conformal invariance follows 
from the conformal invariance of the Weyl tensor 5 , and, second, the 
equation of geometrodynamics, proposed in 1984 in ||, see also [|7| and [§]. 
As the details of the latter are not so well-known, we sum up their properties 
in the Appendix. 

It is the purpose of the present paper to compare these two latter ap- 
proaches. 

2 Notation 

We apply the following notation: Riemann tensor: 

\ _ f a 

aa/3 — J- (3a. 



R X aaf3 - fa - . . . (1) 



Ricci tensor: 

Ra/3 = R a aaf3 ■ (2) 

This is the rule for rearrangement of order of covariant differentiation: 

Y X ;a;(3 = ^ X ;/3;a ~ Y^R^aajSi ^A;a;/3 = ^A;/3;a + Y a R° ' \ a p . (3) 

The Weyl tensor is: 

C a (3^ v = R a /3f_iu — \ [g a ^R/3u + 9/3uRafj, ~ 9auR/3^ ~ 9/3fiRau] + ^ 
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3 The Bach equation 



Bach equation describes the vanishing of the variational derivative of the 
square of the Weyl tensor with respect to the metric. We use it here as given 
in i, i.e.: 

B a ? = B§ + B<$ = . (5) 

Here: 

i 2 1 

B a f3 {R-'i9-) = ~Raf3 ,U ;u + Ra;P;u + Rj3;a;v ~ ^R;»;P + -g a pR ,L ' ;v , (6) 

2 11 

B a /3 {R~> 9-) — ^RRa/3 — ^RauR U P — -^R 2 9aP + ^apR^R^ ■ (7) 

From eqs. (|6|), (|7|) it becomes clear that the quantities B^, depend 
on the Ricci tensor and on the metric. It should be mentioned, that up to 
a multiplication of the whole Bach tensor by a non- vanishing constant, all 
these 9 coefficients in front of the 9 different geometric terms in the right 
hand sides of eqs. @ and (0) are uniquely determined by the conformal and 
covariant properties of the Weyl tensor. 



4 The equation of geometro dynamics 

By the "equation of geometrodynamics" we denote the equation deduced in 
[[]]; it reads (cf. the Appendix for a deduction): 

R a p — -^9apR = —2A a Aj3 — g a pA 2 — 2g a pA l ' ; „ + A a .^ + A^ a , (8) 

where 

A 2 = A^A u g^. 
From the equation (§) we get after taking the divergence 

F a u ]U = , (9) 

where 

F a p = Ap :a — A a $ = Ap. a — A a .p . (10) 
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The trace of equation @ yields: 

R = 6A 2 + 6A». U . (11) 
Then the equation (|8|) can also be written as 

R a p = -2A a A p + 2g aP A 2 + g a pA v . v + A a , p + A p . a . (12) 

5 The definition of the problem 

We apply the following conformal transformation: 

g a p(x)^g' a(3 (x)=g a(3 (x)-e 2 ^K (13) 

The equation of geometrodynamics @ is invariant with respect to conformal 
transformations (|T3|) , if the vector A a (x) is simultaneously transformed as: 

A a (x) -> A' a (x) = A Q (x) - a Q (x) . (14) 
The Bach equation ([5]) has the following solutions: 

— Riemannian spaces with zero Ricci tensor, i.e. with R a p = 0. After a 
conformal transformation we get 

Here (x) is an arbitrary function of the coordinates. Due to the conformal 
invariance of the Bach equation, Riemannian spaces with Ricci tensor (|15|) 
also satisfy the Bach equation. [] 

— Einstein spaces, 



R a f3 = A • g a p, A ^ 0, A = const. , (16) 



1 Hint for a proof: Insert ip — a into eq. (]T^), then eq. ( [l5| ) implies R' a a = 0. The last 
two terms of eq. (^|) coincide, of course, but we wrote them this way we did to ease the 
comparison with eq. (12). 
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as well as Riemannian spaces, that differ from the Einstein spaces by the 
conformal transformation (pT3|). If this a is a constant function a = o"o = 



const., then we get again an Einstein space, this time with A' = Ae 2cr °.0 

As a result, we have two types of conformally invariant equations that 
describe Riemannian space dynamics. Each type of equations can be treated 
as a basis for a conformally invariant theory of gravitation. These theories 
are not equivalent, since there are the following differences between Bach 
equations and equations of geometrodynamics: 

1. The equation (||) is of the 4 th order, while (§) is a second-order equation. 

2. The equation (H) includes the vector A a (x), and the eq. (|^) does not. 

3. For the equation (|8|) a correctly-defined Cauchy problem exists (see J7J). 
As for the equation (|5|), correctness of the Cauchy problem definition has not 
been proved yet.f\ 

4. The equation (|5|) can be derived from the variation principle; no La- 
grangian with the dimension of [length 2 ] exists for the equation (|8|). 

Remember that the Lagrangian for the Bach equation reads: 

L = C aP(a/ C a ^. (17) 
Taking into account the fact that 

LqB = Ra/3fi,uR al3fJ ' l/ — AR a pR a ^ + R? (18) 



2 Solutions of the Bach equation which are not conformally related to an Einstein space 

are called "non-trivial solutions" , there existence is shown in Q , || and Q . 

3 More concretely, there are at least two open questions in this respect: First, that 
point in the usual Cauchy problem formulation of the Bach equation, where one breaks 
the conformal invariance by requiring the curvature scalar to be a constant: it is not clear 
whether this restricts globally the set of solutions or not; second, even for the simplest 
anisotropic cosmological models, namely Bianchi type I, it is not clear whether the state- 
ment "Without loss of generality, one may assume the Bianchi type I metric in diagonal 
form to get the full set of Bianchi type I vacuum solutions." proven for Einstein's theory, 
is valid for the Bach equation, too. This second question shall be discussed more detailed 
in a further paper. 
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is the complete Gauss-Bonnet divergence, the Lagrangian for the equation 
(H) can also be written as Q: 

L = 2R a pR a P — ~i? 2 . (19) 

It is clear that any solution of the Bach equation in the form fll5|) with a 
constant A = (] will also satisfy the equation (||), since (p~5| ) can be reduced 
to (|S|) by equating 

A a (x) = -<p, a (x) . (20) 

The subject of this paper is the solution of the inverse problem, i.e. to find 
out under what conditions solutions of the equation of geometrodynamics 
(H) are also solutions to the Bach equation @. 

6 Solution of the problem 

The Ricci tensor (|l^) will be referred to as geometrodynamical tensor and 
denoted as -R GD . We will show that for R GD there exists the following rela- 
tion: 

B aP (R GD ; g.) = 2F oa ,F p v - l -g a pF^ . (21) 

The proof is lengthy. To make the procedure convenient for check, we will 
divide the proof into 7 steps. 
Step 1. 

In this step R GD from ( |12"1) is substituted to B^l (R..;g..). This yields: 

B§ (R GD -) = ^ (A'") + $ a/3 (AA") + (A' A') . (22) 

The terms' dependence on the vector A a and its derivatives A a ^ is shown 

4 up to surface terms in the action which do not influence the classical field equation 
5 The equation (0) can also be written with a non-zero A-term, and with the trans- 
formations ([[3]), (|l4|), A will be transformed as A — * A' = A • e~ 2a . For our proof it is, 
however, sufficient to assume A = 0. 
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symbolically in brackets. In the explicit form, these terms read: 
&af3 0*4-') — ~ A u .„. a] p — A u . u .p. a + 

J rA a ,v '.p. v + Ap' u . a . u + A v . a .p. u + A u .p. a . u — A a .f3 ,,y . u — Af3- a ,u . v 

$ Q/3 (AA") = 2 0V%) A/3 + 2 (A^, v ) A a - 2g Q(3 A° (A a ' v iV ) - 
—2A a A u p. ll — 2ApA v a . v — 2A"A a .p. v — 2A U Ap. a . u 

(A' A') = iA a;v Ap' v - 2A a , v A v ,p - 2Ap, v A v , a - 
-2A a .p (A". v ) - 2A p . a (A v . v ) - 2g a p (^M w „) 

Step 2. 

In this step, the term $ Q( g (A'") is first transformed with the Bianchi 
identity to the following form: 

(A'") = 2A a B^. ff + 2A", a R^ + 2A v ,pR^, (26) 

then, after the expression for R^ D has been substituted, it reads: 

$ a/3 (A'") = e a p (aa") + e a p (A' a') + e a p (a 2 a') , (27) 

where 



B a p (AA") = 2g af3 A u A° , a , v + 2A V A a ,p, v + 2A v Ap, a , v , (28) 

Q a p (A' A') = 2A a .p (A u . v ) + 2Ap. a (A v . u ) + AA v , a Ar^+ 
+2A v;a Ap> u + 2A v . p A a ' v 

Q a p (A 2 A') = —AA V A a . u Ap - AA u Ap. v A a + 8g a p (MM,.,) - 
-AA v A v . a A p - AA V A V . P + AA a ,pA 2 + AA p . a A 2 

Step 3. 

In this step, all terms are substituted to the relation ( p2| ) for 



(29) 
(30) 



here, (^4 W ) being written in the form fl27|) , where the relations for 



Q a p(AA") ; 9^(A'A0 ; 
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and Q a p (A 2 A') are given by (|28|), (|29|) , (|3CiD , respectively. This yields: 

fl$ (i? GD ; g.) = * a p (AA") + * a/3 (A 2 A') + * a/3 {A A 1 ) , (31) 



where 



q af3 (AA") = 2 (A a ' u - u ) A/3 + 2 (A^. u ) A a - 2g a pA° (A a > v . u ) 
-2A a A v , p , v - 2ApA v ,^ + 2g txP A a A v . v . <t 



V a p (A 2 A') = AA a .pA 2 + AAp. a A 2 + 8g a p (4M"AJ - 
-4^^^^ - 4A"A Al/ A» - ^A u A„. a A p - iA»A v . >/3 A a 



, (32) 



(33) 



V a p {A' A') = AA a , u A p - v + 2A ir>a A v . fj - 2g af3 (A^A^) . (34) 

In this step, the terms with the third-order derivatives disappear from the 
relation for B^p (-R GD ; g ^j ■ The next step will be to exclude the second-order 
derivatives from B^p (-R GD ; g.^j- 

Step 4. 

We transform the expression ty a p (AA"), using the properties of the Rie- 
mann tensor. This yields: 

V a p (AA") = -2g a pA 2 (A v . v ) - Ag a p (AM"AJ - 
—2A a Fp v . v — 2ApFa \v + ZgapA F a v 



If the Ricci tensor is written as (|T2|), then the vector A a shall necessarily 
satisfy the equation (^). Taking this into account, fl35|) will change to 

* Q/3 (AA") = -2^A 2 (.4%) - 4g a p (A X A V A KV ) . (36) 

Note that at this step the terms with second-order derivatives have been 
replaced by the terms with first-order derivative. 
Step 5. 

We will find the resulting relation for B^l 

(^ GD -) = V a p (A 2 A') + Q a p (A' A') , (37) 
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where 

n aP {A 2 A 1 ) = AA a .pA 2 + 4A/3- a A 2 + Ag a/3 (A»A»A^) - 

-AA v A a , v A p - AA v Ap. u A a - AA v A y . a Ap - AA v A v .pA a - 2g aP A 2 {A v , v ) ' 

(38) 

n af3 (A' A') = AA a . v Af3' u + AA v , a A v . p - 2g a0 (A* v A^) . (39) 
Step 6. 

We substitute the expression for R GD in the relation for B^p (-R GD ;g.). 
All algebraic terms (the terms that do not include derivative of the vector 
A a ) are eliminated, only terms with the derivatives of not higher than the 
first order remaining. 

2*3 (R GD ; g.) = ~ Q/3 (A 2 A') + ~ Q/3 (A' A') , (40) 

here 

E af3 (A 2 A') = -Q aP (A 2 A') , (41) 

E a/3 (A' A') = -2A a . u Ap>» - 2A a , v A v , p - 2A p . v A v . a 
-2A v , a A v ,p + g aP {A^A KV ) + g a(i 

Step 7. 

The proof of the relation ([H]). To prove this relation, we must substitute 
the relation @7|) for (i? GD ;#..) and the relation (gg) for (R GD ;g..) 
into the RHS of (^TJ). So, the relation (|2l| ) has been proven. The relation 
(j21~D leads to the following three statements. 

Statement 1 : 

If A a is a gradient vector, i.e. it is in the form (ffij), then the Riemannian 
space with Ricci tensor ( 113 ) is a solution to the equation (^j), since in this 
case F a( g = 0. 

Statement 2 : 

If we take a Lagrangian in the formQ (43) instead of (|T9|), 

L = L + 2F iiV F iXV = 2R aP R al3 - -R 2 + 2F 2 , (43) 

3 



6 This Lagrangian will be referred to as Weyl-Maxwell Lagrangian (WM Lagrangian), 
see also H for it. 
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then the generalized Bach equation will read 



B 'a8 = B al + B 'aj3 = °> 



(44) 




(45) 



The Ricci tensor (|T^) satisfies the equation (|33|). 
Statement 3 : 

If the Lagrangian ( |l9|) is replaced by the Lagrangian 




where const ^ +2, then the Ricci tensor ( |T2"D will only satisfy the dynamic 



equation if the vector A a is in the gradient form. 

7 The spherically symmetric solutions of the 
Bach equation and the geometrodynamics 
equation 

Any spherically symmetric solution of the Bach equation is equivalent (almost 
everywhere and up to a conformal factor) to the Schwarzschild static solution 
in general relativity (GR). Each of these solutions in GR is governed by a 
single parameter, i.e. the gravitational radius. 

For the geometrodynamics equation (^j) a spherically symmetric solution 
can also be given in the static form. However, the spherically symmetric 
solution is now governed by two parameters, rather than one. We will write 
explicitly one of the possible static forms of the spherically symmetric solu- 
tions [] of the equation (^). The peculiarity of this form is that the radial 

7 A detailed description of how to derive different forms of spherically symmetric solu- 
tions to geometrodynamic equation will be covered in a separate publication. 
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component is also the one that describes brightness. 



dr 2 



1+ ir )[}-w^-^) (47) 



+r 2 



d6 2 + sin 2 • d<p 2 



At ro < r << Ro the metric ([17]) is similar to the Schwarzschild metric, 
where ro works as a gravitational radius. It holds: For Rq = 0, metric ([17]) 



is the exact Schwarzschild solution, for ro = 0, metric (|47|) is the exact de 
Sitter solution, but for r ■ R ^ 0, it does not give the Schwarzschild-de 
Sitter solution. 

Meanwhile, the topology of the solution ([17]) is qualitatively different 
from the topology of the Schwarzschild solution. The difference is, that the 
solution ( f4"T|) has a horizon not only at the surface 

r = r , (48) 

but also a singularity of the conformal curvature invariant at the surface 

r = R . (49) 



8 Discussion 

We must point out some facts given by the obtained results. 

• This paper illustrated, with a static spherically symmetric solution as an 
example, the correspondence between the solution of the Bach equation (the 
solution with the vector A a given by A a = —if, a ) and that of the geometro- 
dynamics equation (i.e. the solution with non-gradient vector A a ). The fact 
is that a solution involving a non-gradient vector A a (A a ^ —f,a) cannot 
coincide, in principle, with the Schwarzschild solution in the entire range of 
variation of variables, though in some interval of the radial coordinate it may 
come close to this solution. 

• Solutions of eq. (|8[) describe dissipative phenomena like heat conduction 
and viscosity (see ||). One of the results of this paper is the proven fact that 
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solutions of this type are also solutions to the equation (fP|), that can be 
derived from the holonomic variation principle by a conventional procedure. 
This fact itself is non-trivial. Another amazing feature is that the structure 
of the corresponding holonomic Lagrangian is unambiguous. 

• It is likely that the structure of the Weyl-Maxwell Lagrangian in the 
form ([f3|) enables the construction of an analogue of a Hamiltonian formalism 
for solutions of the equation (g) and to quantize, despite the fact that, this 
equation describes dissipative phenomena. 

• The addition of the Maxwellian term (F^F^) to the Lagrangian (fl3|) 
does not immediately allow to interpret the vector A a as an electromagnetic 
vector potential. The Lagrangian (^) does not yield the Einstein field equa- 
tions with the Maxwell energy-momentum tensor, but the equation (|44l) . The 
Lagrangian (|43|) includes two conformally invariant scalar values with Weyl 
weights, inversely proportional to yj— g. These scalars are: C a f3^uC al3fJ ' l/ and 
F^F^ . They give a complete list of the Weyl space invariants 0. In fact, 
the Lagrangian that includes both invariants is the only possible one among 
those able to produce a dimensionless action with zero Weyl weight, in the 
case of Lagrangian being constructed immediately for the Weyl space. For 
this reason this represents a more natural interpretation for the vector A a as 
being a Weyl vector. 

Appendix: Deduction of the properties of ge- 
ometrodynamics 

Here we give a deduction of basic facts connected with the equation of ge- 

ometrodynamics. This deduction is both a new approach and also more 

general than that ones given in the literature up to now. 

Assume that the source for the Einstein field equation shall be composed 

8 The quantities E a(}llv C aSia ''C a p llv and E^^F " 3 ' F» u have the same Weyl weight, 
they are, however, pseudoscalars. Here, E.... denotes the dual of the Weyl tensor. 
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from a vector field A a , where at most first derivatives of it shall enter the 
energy-momentum tensor. A possible A-term can be omitted without loss of 
generality, | and terms of cubic and higher degree in A a shall be omitted. 
Then we are left with the following equation 
R 

Ra/3 - -^g a /3 = ci A a A(i + c 2 A(a-p) + c 3 g a pA 2 + c 4 g a pA^. a (50) 

containing the four free parameters c$, i = 1, ... ,4. Here, A 2 = A 1 A 1 and 
the round symmetrization bracket is defined via 

A {a .p } = - {A a .p + Ap. a ) . (51) 

This equation (|50| ) goes over to eq. (||), if the quadruple q takes the form 

d = (-2, 2, -1, -2). (52) 

The trace of eq. flSOD reads 

-R=(ci + Ac 3 )A 2 + (c 2 + 4c 4 )A 7 . 7 . (53) 

If one inserts eq. ( |52| ) into eq. (|53|), one gets, of course, eq. (|TT|). 

We take the covariant divergence of eq. (^0]). By use of the Bianchi 
identity, the l.h.s. identically vanishes, so we get 

= ciAyLp + c x A^A p .^ + + ^UA P + 2c^A r3 + c A A] lP . (54) 

For the further calculation it proves useful to define a second vector field, Tp, 
via the equation 

Tp = A a R a p . (55) 
With eqs. (0) and Q53|) we find 

Tp = c 2 A<*A ia .,p } + Ap ((| - c 3 )A 2 - (| + c±)A\) ■ (56) 
To keep simplicity, we again require the absence of cubic terms, i.e. 

ci = 2c 3 . (57) 

9 With A 7^ the result would be analogous. 
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For the quadruple eq. ( (52]) this condition is fulfilled. 
Now we apply eq. (|3|) and find with eq. ([55]) 



A a 0a - A a ]ap = Tp . (58) 

Eq. (|56p with eq. ( [57D gives 

Tp = c 2 A a A {a;P} - (| + Ci )A p A\ n . (59) 

Eq. |>|) with eq. (0) gives 

= 2c 3 AlA + M 7 ^) + + f n ^/3 + c 4 ^ 7/3 . (60) 

The goal is now to find out under which circumstances eq. (^) can be deduced 
from eqs. Q58|) till Eq. (|) with eq. (|10D reads 

^A P = A\. (61) 



Therefore, the term A y * of eq. (|60T ) is superfluous. We take the three 
eqs. (|58|) till (|60|), cancel the vector Tp and the term A', Lg and get afterwards: 



= 2c 3 ^ 7 ^ + 4c 3 AM (/3;7) + |A 7 ;/37 + ^UA P + 

+c 4 A! iPl - czc^A^ + c 4 (| + c 4 ) . (62) 



To get a relation like eq. (|61|), the nonlinear terms in eq. ( |62"D have to be 
cancelled. This takes place if 

4c 3 = c 2 • c 4 (63) 

and 

2c 3 = -C4 (| + cA . (64) 

C2 = implies q = for all i, which is a non-interesting case. Therefore, 
we assume c 2 7^ 0. We are free to redefine the vector Ap by multiplying 
it with any non- vanishing constant, because we can compensate this by an 
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appropriate redefinition of the constants q, see eq. (|50|). We use this freedom 
to choose the value C2 from eq. fl52"|), i.e., C2 = 2. Then eqs. fl63]), ( |64[ ) read: 



4c 3 = 2c 4 ; 2c 3 = -c 4 (l + c 4 ) . 
The solution with c 4 = 0, which implies also c\ = C3 = 0, turns eq. ( |60"D into 

This special case shall now be omitted, then we get as the only other solution 
c 4 = —2, c 3 = —1, and with eq. (|57]) finally Ci = —2. This is together just 



the original quadruple eq. (|52|). So, we have shown, that eq. (£]]) indeed 
follows without further assumptions from the divergence of eq. (§), and we 
have shown in which sense the coefficients q are uniquely determined. 
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